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1 . INTRODUCTION 

From the sensitivity studies performed by Baker and Paegle (1983) and 
Paegle and Baker (1983) with the Goddard Laboratory for Atmospheres (GLA) 
analysis/forecast system (Kalnay et al . . 1983), it was revealed that the 

forecast errors in the tropics affect the forecastability of midlatitudes in 
some cases. Apparently, the forecast errors occurring in the tropics can 
propagate to midlatitudes. Therefore, the systematic error analysis of the 
GLA forecast system becomes a necessary step in improving the model's forecast 
performance. Our major effort of this study is to examine the possible 
impact of the hydrological-cycle forecast error on dynamical fields in the 
GLA forecast system. 



2. FORECAST EXPERIMENTS 


Twenty-eight 5-day forecasts were performed every two days for both 
Special Observating Periods (SOP). The forecast experiments of S0P1 
started from OOOOGMT January 5, 1979 and those of S0P2 from OOOOGMT May 1, 
1979. The initial conditions of these forecast experiments are provided by 
the GLA FGGE Ill-b analyses (Baker, 1983). The analysis scheme employs a 
successive correction method of the Cressman (1959) type, modified to account 
for differences in data density and statistical estimates of the error 
structure of the observations. Nonlinear normal mode initialization is not 
used in the GLA forecast system, but a Matsuno (1966) time-differencing 
scheme is adopted to control the inertial gravity waves excited by the 
imbalance between the mass and motion field. 


The ensemble averages of 28 forecasts will be used to illustrate the 
forecast errors. 

3. SURFACE BALANCE AND HYDROLOGICAL CYCLE 

After Day-1, the surface heat budget of forecasts (Table la) becomes 
relatively stable in time for both SOP’s. Although the surface shortwave and 
longwave radiation are larger than their annual-mean climatic values, the sum 
of these two surface radiation components results in a smaller net surface 
radiation than climate. Climatologically, the surface (sensible + latent) 
heat flux should be balanced by the net surface radiation. By taking 
computational errors into account, the model’s net surface heat flux is 
small . 

The GLA hydrological cycle forecasts of both SOP’s are shown in Table 


lb. The global-mean evaporation over the entire forecast period is slightly 
smaller than the climate inferred from Jaeger’s (1983) estimation of 
precipitation. The precipitation forecasts significantly exceed the 
evaporation forecasts. Even the hydrological cycle of Day-5 forecasts still 
does not reach equilibrium. During the spinup of the hydrological-cycle 
forecasts, the S0P2 precipitation forecasts exhibit an overshooting at Day-1. 
The GLA forecast system does not employ the nonlinear normal mode 
initialization nor the shallow cumulus convection. The Arakawa’s (1972) 
cumulus convection scheme is used in this model. However, the characteristics 
of the GLA hydrological cycle forecasts are also common to the forecast model 
of the European Center for Medium Range Weather Forecast (ECMWF) analyzed by 
Heckley ( 1985 ). 

4. FORECAST ERRORS OF ZONAL-MEAN STATES 

Since the forecasted precipitation exceeds the forecasted evaporation, 

it is expected that the forecasts would become too dry. As shown by the 

Day-5 forecast errors of the zonal-mean specific humidity (q ) (Figure la) 

z 

during SOP2, the middle and lower troposphere are dried out and the planetary 

boundary layer is moistened. The major precipitation exists in the tropics. 

The overforecasted precipitation of both SOP’s also occurs there. Therefore, 

a warm bias of the zonal-mean temperature (T ) forecasts may be caused by the 

z 

latent heat released by the overforecasted cumulus convection activities over 

the tropics. The Day-3 T forecast errors during SOP2 (Figure lb) confirms 

z 

our argument. The north-south T gradients are also overforecasted, 

z 

particularly’ the midlatitudes of the winter hemisphere. According to the 

thermal wind relationships, the zonal-mean zonal wind (U ) should also be 

z 

overforecasted in midlatitudes. As exhibited in Figure lc, the Day-3 U z 

forecast errors during SOP2 behave in a relatively systematic manner. The 
low-level tropical easterlies and midlatitude westerlies are overforecasted, 
while the high-latitude westerlies are underforecasted. 

Our previous study (Chen and Baker, 1986) suggested that the atmospheric 
diabatic heating maintains the tropical divergent circulation. The time 
evolution of the forecasted tropical divergent circulation may follow that of 
the precipitation forecasts shown in Table 1. As revealed from the Hadley 

cell portrayed by the mass flux function = f 1000 m ' 5 v dp) forecasts 

z z 

p 

(Figure 2), the overshooting of precipitation forecasts reflects in the 

maximum intensity of Hadley cell at Day-1 forecasts. By comparing Day-1 and 

Day-5 forecasts, the intensity of Hadley cell becomes weaker when the forecast 

time extends. To be more quantitative, the forecast errors of the zonal-mean 

meridional wind (V ) show that the southerly and northerly anomalies appear, 

z 

respectively, in the lower and upper tropical troposphere of Day-1 forecasts. 
Evidently, Hadley cell is overforecasted at Day-1. In contrast, the 
upper-level northerly anomalies are replaced by southerly and the lower-level 
southerly anomalies disappear in Day-5 forecasts. This implies that the upper 
branch of Hadley cell is underforecasted at this time. 


Except for the precipitation and $ overshooting in the SOP2 forecasts, 

z 

the salient features of the q_, u , T and V„ forecast errors are also common 

z z z z 

to SOP1 forecasts. 

5. FORECAST ERRORS OF HORIZONTAL FIELDS 

In accordance with the zonal-mean divergent circulation depicted above, 
the 200-mb divergent circulation [(Vp.Xy) (200 mb)] forecasts of SOP2 are 

shown in Figure 3. At Day-1 forecasts, three divergence centers are located 
over three tropical continents. The most dominant is the divergent 
circulation associated with the monsoon region. The clear divergent center 
over Central and northern South America at Day-1 forecasts is replaced by a 
relatively weak convergence center at Day-5 forecasts. This significant 
change of divergent circulation explains why the upper branch of Hadley cell 
is underforecasted at this time. Some previous studies, e.g., Silva Dias et 

al , (1983) demonstrate that the diurnal cycle of cumulus convection is vital 

in maintaining Bolivian high. It seems that the diurnal cycle of the GLA 
cumulus convection scheme over this region does not effectively play its 
designated role. The [(Vp,Xy)(200 mb)] forecast errors during SOP1 exhibit 

the underforecast of east-west Valker circulation. 

The major precipitable water (W) exists in low latitudes, particularly 
over three tropical continents and along the ITCZ. The Day-5 W forecast 
errors (Figure 4a) show that the W underforecast previously pointed out occurs 

over the entire tropical global belt between about 30°S and 30°N. The 

significant W forecast error over high mountains may be attributed to 
erroneous data. The W underforecast is caused by the precipitation 
overforecast. The water balance equation 

<V.Q d > = <V 2 x q > = E - P (1) 

indicates that the precipitation (P) overforecasts should be supported by the 
overforecast of convergence of water vapor flux. To facilitate discussion, 
the divergent component (Q^) and potential function (Xq) of water vapor flux 

are introduced in Eq. (1). The Day-1 forecast errors of (Qq.Xq) (not shown) 

show that the precipitation overshooting corresponds to the overforecast of 
convergence of water vapor flux toward tropical continents and ITCZ where the 
major precipitation exists. After this overshooting, the convergence of water 
vapor flux toward the Asian monsoon region and the region around Central and 
northern South America is underforecasted (Figure 4d). This drastic (Qq.Xq) 

forecast change agrees with that of (Vp,Xy)(200 mb) shown previously. During 

SOP1, the W forecast errors behave similarly, as those during SOP2 and the 
( vv forecast errors correspond to the underforecast of east-west Walker 

circulation. 


Yanai et al. (1972) demonstrated that the tropical water vapor and heat 


budgets are closely related to each other through the latent heat. The 
positive values of T(300 mb) forecast error (Figure 4c) and tropical 

convergences of (Qq.Xq) forecast error for S0P2 Day-5 forecasts are 

relatively consistent. This consistency indicates that the excessive latent 
heat released due to the precipitation overforecasts is responsible for these 
T(300 mb) forecast errors. 

The comparison between the SOP2 Day-5 U(200 mb) forecasts (not shown) 
and the forecast errors (Figure 4d) shows that either jet streams themselves 
or the equatorward sides of jet streams are overforecasted. The contrast 
between positive forecast errors of T(300 mb) and U(200 mb) (Fig. 4d) that the 
thermal wind relationship holds between them indicates this. 

During the forecast spinup, the thermal wind relationship is reached 
through the adjustment by ageostrophic circulation. This physical adjustment 
can be realized through the generation of kinetic energy by ageostrophic flow 
(-V»V<£) which is roughly balanced by the divergence of kinetic energy flux 
[-V»(Vk)l. Their relationship is essentially the kinetic energy equation 

0~-V.(Vk) D -V.Vtf> = -V 2 *^ - vTV? . (2) 

The divergent component t(Vk)p] and potential function (Xy^) of kinetic 

energy flux are used. The kinetic energy (k) is generated in the upstream 
side and diverged to the downstream side of the jet where the kinetic energy 
is destroyed by ageostrophic effect. The SOP1 Day-5 forecast errors of k(200 
mb) (Figure 5a) are consistent with the observation of the U(200 mb) forecast 
errors shown in Figure 4c. The SOP1 Day-5 forecast errors of (Vk)p(200 mb) 

(Figure 5b) are opposite to the SOP1 Day-5 (Vk)p(200 mb) forecasts (not 

shown). This situation indicates that the thermal field forecast errors 
possibly alter the ageostrophic circulation forecasts in such a way that the 
generation of kinetic energy in the upstream side is underforecasted. 
Therefore, there is less kinetic energy transported to the downstream sides of 
jet streams to be destroyed. 

The tropical divergent circulation is maintained by the diabatic 
heating. In turn, the precipitation forecast errors can cause the divergent 
circulation forecast errors. The planetary-scale interaction between 
divergent (Vp) and rotational (V R ) may offer another possible adjustment 

process causing the forecast errors of jet streams. Following Chen et al . 

(1988), the interaction between these two flow components can be understood 
through the kinetic energy equations of divergent (kp) and rotational (k R ) 

flows: 


0 ~ - (f+£Ku R Vp-UpV R ) - Vp*V<£ , (kp equation) 
C(kD,k R ) 


( 3 ) 


(k R equation) 


(4) 


0 " - V*(V R k) + (f+O^Vp-UjjVj^) - V R .ty 

C(kjj,k R ) 


The SOP2 Day-5 forecast errors of Cdcp,^) are shown in Figure 6. The 

,k R )(200 mb) forecast mean that the excessive kp(200 

mb) is converted to maintain the excessive k R (200 mb). Some consistency 

,k R )(200 mb) forecast errors indicates 

that the U(200 mb) forecast errors are possibly caused by those of divergent 
circulations. 


existing between the U(200 mb) and C(kp 


positive errors of CKk^ 


5. CONCLUDING REMARKS 

The GLA forecast system does not employ nonlinear normal mode 
initialization nor the shallow cumulus convection scheme. The cumulus 
convection scheme proposed by Arakawa (1972) is adopted in this model. The 
systematic errors of thermal and dynamical fields in the GLA forecast system 
differ from those of the ECMWF forecast model, but the GLA’s deficiencies of 
moisture and precipitation forecast are common in the ECMWF forecast model 
(Heck ley, 1985). According to Tiedtke et al . (1988), the incorporation of a 

shallow cumulus convection is most vital to the improvement of the ECMWF 
forecast performance over the tropics. The impact of the shallow cumulus 
convection on forecasts is one of the future efforts which should be pursued 
in the GLA forecast system. 
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(b) cpz forecast (SOP2~Day5) 


(c) Uz forecast errors (SOP2-Day3) 


(c) Vz forecast errors (SOP2-Day1) 


(d) Vz forecast errors (SOP2-Day5) 


The zonal-mean cross-sections of 

100 

S0P2 Day-3 forecast errors 

of 


a) specific humidity (q ), 


150 

b) temperature (T ) and c) 
z 

zonal 

200 

250 

wind (U z ) . 


300 

400 

500 

ORIGINAL PAGE IS 


700 

OF POOR QUALITY 


850 




A \ 


60S 30S E 


n h ■ \ r-XS 


30N 60N 9C 


a) except for S0P2 Day-5 forecast, c) S0P2 Day-1 forecast errors 

of zonal-mean meridional wind (V ) , and d) same as c) except for 
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S0P2 Day-5 forecast. 






































